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abstract A mean free path of nucleon moving through nuclear 
matter with kinetic energy of more than lOOMeV is formulated based on the 
bare vertex nuclear Schwinger-Dyson (BNSD) method in the Walecka model. 

The self-energy which is derived from the higher order diagrams more than the 
forth order includes the Feynman part of propagator of energetic nucleon and 
grows up rapidly as an increase of kinetic energy. To avoid too large growth 
of these diagrams, meson propagators are modified by introducing some form 
factors to take account of a internal structure of hadron. It is confirmed that the 
mean free path calculated by the BNSD method agrees good with experimental 
data if a reasonable form factor is chosen, i.e., a dipole (quadrupole) type of 
form factor with a cut-off parameter about 750 MeV ~ 1000 MeV (1200 MeV 
~ 1500 MeV). 
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1. INTRODUCTION 


The nucleon mean free path in nuclear medium is one of the important 
quantity to show the property of nuclear matter as well as the binding energy 
and the incompressibility etc. It is experimentally extracted from the proton (or 
neutron)-nucleus total cross sections for many targets at a spread of energy [1] [2]. 

Recently there have been the numerical results of the nuclear mean free path 
in nuclear matter in the framework of relativistic approach based on the Dirac 
phenomenology [3] , in the relativistic impulse approximation [4] [5] and in the 
density dependent Hartree-Fock approximation [6].In this paper we evaluate the 
nucleon mean free path in nuclear matter in the a- u> model based on the nuclear 
Schwinger-Dyson formalism . 

In the last several years ,the two loop correction and the ring-sum correction 
were calculated in the framework of renormalizable a — w model and these 
contributions to the energy density are much in the Feynman part and in the 
finite density part [7] [8] [9] [10]. It is known that the large contributions of the 
Feynman part are reduced by introducing form factors at all vertex in the loop 
diagram [11] [12] or by vertex corrections[13]. These procedures for reducing 
short-range distance contributions are qualitatively acceptable at present to 
take account of the internal structure of hadron. By using the bare vertex 
nuclear Schwinger-Dyson (BNSD) method we showed that the loop corrections 
in the density part were also hardly reduced because the these corrections were 
strongly canceled among the contributions of three components ,i.e., a- and 
w-mesons and a-tu mixture [14] [15]. 

Until now we have studied about the saturation property of ground state 
of nuclear matter based on the BNSD method[14][15]. In these studies, we 
determined two meson-nucleon coupling constants to get the minimum energy 
at the normal density, taking care to remove the instability near the normal 
density by using some recipes [16] [17]. After this work, we studied optical 
potentials for nucleon near Fermi surface under the BNSD approximation by 
using the results of nuclear matter[17], and we obtained the reasonable optical 
potential both in the real and imaginary parts compared with empirical findings 
[18] [19].. 

In this paper we estimate the mean free path of nucleon which travels through 
nuclear matter having the kinetic energy of more than 100 MeV above the Fermi 
surface. We will show that an optical potential for such a high energy nucleon 
outside the Fermi sphere is derived from the Feynman part of self-energy which 
corresponds to the higher order diagrams more than the fourth-order, besides 
the density part of self-energy. These higher order diagrams yield very large 
contributions to the optical potential both in the real and imaginary parts as well 
as large contributions of the Feynman part of energy density of nuclear matter. 
We point out that there need meson propagators modified by introducing some 
form factors to make the mean free path calculated by the BNSD method agree 
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with ones extracted in this kinetic energy region. 

The organization of this paper is as follows. We develop the BNSD method 
for a energetic nucleon traveling far above the Fermi surface in Sec. II. The 
numerical results and discussion are shown in Sec. III. We summarize our work 
of this paper in Sec.IV. 

2. FORMULATION 

We adopt the Walecka model[20] which consists of three fields, the nucleom/),the 
scalar cr-meson <f> and the vector w-meson V^.The lagrangian density is given by 

L = -^( 7 ^ + M)ip - i(<9 + m 2 s (j) 2 ) 

+ ^mlVuVft) + gs^ip(j> + ig v '4>lv,Wn, (!) 

where = d^V u — <9„V^ and M, m s , m v , g s and g v are nucleon mass, cr-meson 
mass, w-meson mass, cr-nucleon and w-nucleon coupling-constants respectively. 
The nucleon propagator is obtained by following form, 

G{k) = G F {k) + G D (k) 


— 1 77r 

- + + M* k )—6(k* 0 - E* k )0{k F - |k|), 


+ M* - ie 


( 2 ) 


E* k = J k * 2 


Mf, 


(3) 


M* = M + E s {k), 


(4) 


K = (k*,ik* 0 ) 


= ^k(l + E v (k)),i(k 0 + E 0 (k))j , (5) 

where G F and G F are the Feynman part and the density-dependent part, re¬ 
spectively. 

The nucleon self-energy £(fc) is classified as shown in Fig.l 

Fig-1 


£(fc) = + E SDD (fc) + E SDJ ?(fc), (6) 
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where £//, £ sdd and £ sdf are the Hartree term ,the SD density term and the 
SD Feynman term respectively. And this nucleon self-energy is alternatively 
classified according to the property of Lorentz transformation, 

£ = £ s (k) - 70 £ 0 (fc) + iji ■ kiT, v (k), (7) 

where £ s ,£ 0 and £.„ mean the scalar component ,the timelike component of 
vector type and the spacelike component of vector type. 

In the BNSD method, the meson propagator satisfies the following Dyson 
equation as shown in Fig.l, 

D{q) = D 0 (q) + D 0 (qMq)D(q), (8) 

where Do(q) denotes the free meson propagator and 11(g) denotes the meson 
self-energy. Then the meson propagator has the real and the imaginary parts, 

D(q) = ReD(q) + Hm.D(q). (9) 

To derive the imaginary part of self-energy we execute a Wick rotation go —» 
igo in the go-integral of the SD Feynman term and rewrite T,sdf as follows, 

XSDF = E ° s n DF + £ SDF, (10) 

ReEgDF = ReY.°s n DF + R^°sd F , ( 11 ) 


ImZsDF = ImT,° s n DF . (12) 

The superscripts ”on” or ’’off” of T,sdf denote that the nucleon propagator 
included in the SD diagram is on-shell or off-shell, respectively. We note that 
the real part of T,sdf is composed of two terms, i.e., a finite term ReT,g n DF and 
an infinite term R eY,°J FF , while the imaginary part of T,sdf is only one term, a 
finite ImEgjjp. We drops the infinite real self-energy as stated in the previous 
paper [17]. 

We obtain the explicit expression of real and imaginary nucleon self-energies 
corresponding to the Feynman diagrams given in Fig.2 as follows, 

Fig-2 

2 g 2 g 2 [ kF „ M* f 1 

i?eE s (fc) =- q ~£ dq / dxA o(R) 

7T 2 mi 8 tH J 0 E* J_ ± 

a 2 r kp M* r 1 

+ s^l "e 

~ 2 pk pi d 2 

+ o2/ q2dq dx-^[^-(A 0 n m D m ) R +(a 0 h s d s ) r ] 

^ Jk F J-l Rq R 2 
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9v 

8n 2 


M* 


R 2 


R 2 


’ k,F 


q 2 dq I dx—^-[-4:(D{)IItDt)R+—{—(DQlliDi)R+(DoHtDt)R— — (DoTl rn D rn )R}} 

£j q 


f-i 


R 2 

rl p 2 


+ ^l/L dx w ReDM 

dq J dxImD s (R) 


R 2 

(13a) 


ImT, s (k ) = 


f k 2 M q 

/ 9 "f*" 

'fc F *4 


87T 2 


'&.F 


M* 

J 2 -?-dq 

E* \ 


dx{4:Im.Di(R) + (3 + 


S%)(ImD t (R) - ImD^R))} 


+ T^ f 1 2 dq [ dxImD m (R )^-, (136) 

4tt 2 J kp i? 2 

O „2 7.3 2 /-fc F /•! 

i?eSo(fc) = — ~2^2~& + / <?dq dxAo(R) 

tt z to 2 3 87t 2 Jo 

+ ^ q2dq j\ dx2D o( R ) 

(.2 pk nl j ^2 

q -~2 / 1 2 dq dx[^-(A 0 U m D m ) R + (A 0 H s D s ) r \ 

^ Jfc F J _1 l ? 2 


J 2 pk />1 ^>2 ^»2 

+ / rfa:[—2(Dontl?t)H—(Z?on;H;)ii+(£>ontr> t )fl—^-(£>on TO £> m )/{}] 

ott Jfc F J-1 R 2 R 2 


^9s9v 

8 tt 2 




M* i ? 2 

Jx—— ^r- ReD m (R ), 
Eg R 2 mV h 


(14a) 


ImT,o(k) 


9 s 
87T 2 




q 2 dq / dxImDs(R) 


i-i 


Q 2 fk /*! d 2 

—^ / q 2 dq / dx{2ImDi(R) + (2 + -=^)(ImD t (R) — ImDi(R))} 
^ Jk F J-i R 2 


g s gv 

47T 2 




M* i? 2 

dx ~EA IrnDm ( R ) : 


(146) 
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n 2 r kp f 1 a* hr 

n 2 f k r 1 a* hr ft 2 

-g% / q 2 dq dx^M A 0 n m D m ) fi + (A 0 n s D 5 ) K ] 

J k jfr 1 / — 1 ^ Jtb 

n 2 i k r 1 n*kr R 2 R 2 

+^TT 2 / [-2(£> 0 n t A)fl-^r{-(i>onzA)fl+(r>on t D t )fl-^-(r>on m i? m )fl}] 

87r Jfc F 7_i i? 2 it 2 

(15a) 

r / 2 /■* r 1 r/*A-r 

/mS » w = - toW., L, d *^r I ™D.(K) 

n 2 r k f 1 n*k r R 2 

-g^a J k q 2 dqJ^dx^{2ImD l (R) + (2-^)(ImD t (R)-ImD l (R))}, 

(156) 

where R = k — q, and p s and ps denote the scalar and the baryon densities 
respectively, and Ao and Z?o denote the free propagators of cr- and w-mesons, 
respectively. The subscripts s,l,t and m of D(i?) and 11(7?) denote the compo¬ 
nent of a-meson, the longitudinal and transverse components of w-meson and 
the component of mixture of cr- and w-mesons, respectively. The detailed ex¬ 
pressions of D(R ) are given in Ref. [16] and the analytical expressions of 11(7?) 
are given in Ref. [10]. 

The last three terms in Eqs.(13a) and (14a) and the last two terms in 
Eq.(15a) are contributed from the same higher order diagrams which yield the 
imaginary part. 

We modify the meson propagators in the SD Feynman part of self-energy 
by introducing some form factors at the vertices in agreement with the point of 
view in Ref. [11] [12] [21] to examine the short range interactions and the finite 
size of hadrons ,as follows 

D(q 2 ) -> [F(q 2 )] 2 ■ D(q 2 ), (16) 

where F{q 2 ) = 1/[1 + g 2 /A 2 ]™ . The case of n=l is a monopole type of form 
factor and n =2 is a dipole type, etc. 

Optical potentials are defined from self-energies E as 

u s = = Usr + iUsi , (17) 

1 + 2-i v 
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Uy = —-, 0 v —-=UvR + iUvi , (18) 

where E is the energy of objective nucleon propagating with the momentum k 
satisfying the following dispersion relation, 


E = Vk 2 + (M + U s ) 2 + U v . 

Eq.(19) is rewritten as 


— + V + iW = E-M + 


(E - M) 2 
2 M 


with Schrodinger equivalent potential form 


(19) 


( 20 ) 


V — Usr + Uvr H —— -l \ r + ~2m^sr + Uv! - Ugj - Uy R ), (21) 

(E - M) 1 

W = Usi + U V i H-—- U VI + —( UrrUsi ~ U VR U VI ). (22) 

The nucleon momentum is complex as well as the optical potentials and can be 
expressed as 


|kj = k.R + iki. (23) 

and then we obtain the nucleon mean free path as follows, 

A= Wi 


1 

2 


-m(e-m-v + ^E 


+M 


(e-m-v+ 


(E - M) 2 
2 M 



(24) 


3. RESULTS AND DISCUSSION 


In this section we calculate self-energies of the nucleon propagating with the 
energy E and the momentum k , convert them into the optical potentials in 
the Schrodinger equivalent form , and evaluate the nucleon mean free path from 
Eq.(24). The relation between the energy E and the momentum k is given by 
approximating Eq.(19) as follows, 

E = y/k 2 + (M + ReT, s (E , k)) 2 - ReT, 0 (E , k). (25) 


1/2 
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The nucleon momentum k, the upper limit of integral in the real and imaginary 
self-energy, is determined from Eq.(25) by putting the nucleon energy E = 
Em + M ,where is the nucleon incident energy. 

In this paper we determine the coupling constants of the cr-nucleon and 
cc-nucleon to satisfy the saturation property of nuclear matter at the normal 
density p=0.170 fm -3 (kp = 1.36 fm _1 ) , and 0.193 /m -3 (kp = 1.42 /m _1 ) 
and summarize the parameters in Table. 

Table 

Using these parameters (we choose the parameters in the case of the normal 
density p= 0.17 fm -3 ), we calculate the nucleon mean free path. We start with 
evaluating the higher order diagrams more than the fourth order in the nucleon 
self-energy shown in Fig.2. Since these diagrams include the Feynman propa¬ 
gator of nucleon on-shell in the intermediate state, we are afraid whether these 
diagrams give a large contribution to the self-energy. In Fig.3, we show contri¬ 
butions of cr-meson, to- meson and a-u; mixture to ReE s and ReEo, respectively. 
As the increase of kinetic energy of nucleon, E — M, each component, E CT , E w 
and S' 7- " , grows almost linearly, because the meson propagator ~ 1 /q 2 and 
the Jacobian q 2 factor cancel each other in q- integrals of self-energies with the 
upper limit of k. In ReE s , additive contributions of cr-meson and a-u> mixture 
cancel strongly with the one of w-meson and as a result the net contribution is 
small. On the other hand, contributions of three components to ReEo are all 
additive and as a result the self-energy is very huge and the net contribution 
becomes nearly equal with the Hartree contribution. 

Fig.3(a), Fig.3(b) 

In Fig.4, we also show the contributions of three components to ImE s and 
IrnSo, respectively. The contributions of cr-meson and a-co mixture to ImEo 
are far smaller than the one of w-meson. The w-meson dominance should be 
marked. The value of ImEo at E — M = 200MeV are considerably reduced 
in comparison with the one derived from the fourth order diagrams in Ref. [21]. 
The higher order contribution more than the fourth order is very important [17]. 
The value of ImE,, is also reduced at the same kinetic energy and, however, has 
the negative sign when E — M > 350 MeV, which is opposite to the sign of ImE s 
obtained from the fourth order diagram. The change of sign is easily understood 
from Fig.4(b). In lower kinetic energy region ImE s is positive because the 
contribution of ui meson is large and positive compared with the small and 
negative contributions of others. On the other hand, in higher kinetic energy 
region, the sign of ImE, changes to the negative one because the contribution of 
u> meson changes to the negative one and becomes large. It should be remarked 
that ImE originates from the Feynman part of self-energy. So, although ImE is 
small if E — M < 100 MeV, it grows up unphysically if E—M > 100 MeV as well 
as the vacuum effect for the energy density in nuclear matter is unphysically 
large. Then, we introduce form factors into each vertex by modifying meson 
propagators, D(q) —> [ f(q 2 )] 2 • D{q) as shown in Ref.[21]. It is noted that, 
since the SD diagrams are composed of full meson propagators with a ring-sum 



correction, a dipole (quadrupole) type of form factor introduced into the SD 
diagram corresponds to a monopole (dipole) type of form factor into the fourth 
order diagram. 

Fig.4(a), Fig.4(b) 

There are two kinds of T,sdf leading to the real self-energy. One is the 
above-discussed self-energy derived from the same diagram which yields the 
imaginary part. The other is the self-energy derived from the diagrams including 
the Feynman propagator of nucleon off-shell. The renormalization procedure 
showed that ReE sdf is very large. So we must develop a new recipe which 
takes account of the size of hadron. Then, we should evaluate the two diagrams 
at the same time when we want to know their contributions to the real self¬ 
energy. In the present situation, therefore, we do not pick up both of them 
though we are afraid of the violation of the dispersion relation between the real 
part and the imaginary part. 

In Fig.5(a) and (b), we show the imaginary part and the real part of Schrodinger 
equivalent potentials , in cases of two types of form factor, the dipole with A 
= 750 (solid curve), 1000 (dotted curve) and 1500 MeV (dashed curve) and 
the quadrupole with A = 1200 (bold solid curve) and 1500 MeV (bold dotted 
curve). The solid dots are empirical information [18][19]. As the kinetic energy 
increases, the imaginary potential decreases slowly if E — M < 300 MeV and 
rapidly if E — M > 350M eV. This feature is particularly in character with the 
dashed curve which is closer to the one without the form factor. The reason is as 
follows. In the rough, the imaginary potential is proportional to ImE s — ImE 0 
and ImE s changes from positive sign to negative sign at E — M ~ 350 MeV. 
The numerical data of the dipole type of form factor with A = 750 MeV almost 
correspond to the data of the quadrupole type with A = 1200 MeV. Similarly, 
the dipole type with A = 1000 MeV corresponds to the quadrupole type with 
A = 1500 MeV. Our numerical results are a little different from empirical infor¬ 
mation in the tendency in the low kinetic energy region. 

Fig.5(a), Fig.5(b) 

The nucleon mean free path is written by using Schrodinger equivalent po¬ 
tential form as given by Eq.(24). In Fig.6 we show the mean free path in the 
cases of some form factors. When we choose the dipole type of form factor with 
cut off A=750 MeV and the quadrupole type of one with A=1200 MeV , the 
values of mean free path are good agreement with the experimental ones, but 
the experimental data are extracted from the reaction of proton and nuclei and 
so we consider that the reasonable form factor is the dipole type with A cs 750 
~ 1000 MeV or the quadrupole type with A ~ 1200 ~ 1500 MeV for nuclear 
matter. 

Fig.6 

Furthermore, we calculated the similar calculations using coupling constants 
determined at the normal density p=0.193 (fcp=1.42 fm -1 ), but the results were 
nearly equal. 
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4. SUMMARY 


We evaluated the nucleon mean free path of the energetic nucleon based on 
the BNSD method and compare it with the ones extracted experimentally [1], 
and took a good agreement in the case of modification of meson propagators 
by introducing dipole type of form factor A = 750 MeV ~ 1000 MeV and 
quadrupole type with A = 1200 MeV ~ 1500 MeV. We confirmed that the 
BNSD method was useful for the derivation of optical potential for a nucleon 
traveling far above the Fermi surface if we consider the finite size of nucleon or 
some effects of the short-range corrections of nuclear force. 

The real self-energy ReEg^ grows unphysically as the increase of kinetic 
energy. Even if it can be reasonably reduced by modifying meson propagators, 
the residual Feynman self-energy ReEg^f F is expected also very large. We assure 
again the conclusion discussed in the previous paper [17] that we do not pick 
up Re'Eg n DF when we drop out ReE °/ FF - On the other hand the imaginary part 
of E sdf is only ImE^, F . The cut off parameter of form factor introduced into 
ImE gjjp looks somewhat small but reasonable when we regard this parameter 
as a parameter of nucleon size. 

The real part of optical potential in the Schrodinger equivalent form increases 
linearly as the increase of kinetic energy, keeping a good agreement with the 
experimental findings in the lower kinetic energy. 

There remains a problem in the imaginary part of optical potential that its 
magnitude is too small in the lower kinetic energy region when a form factor is 
chosen to obtain the experimental data of mean free path in the higher kinetic 
energy region. In our previous work and this work, we took account of the 
vacuum part of meson self-energy with the low 4-momentum transfer to remove 
a instability around the normal density. The vacuum part of meson self-energy 
is also inserted into the denominator of the meson propagator in the integral 
of ImEg^. The vacuum part in this case has the high 4-momentum transfer 
and so is expected too much. As our next work, we will study the role of the 
vacuum effect for meson self-energy in the imaginary potential. 
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Table and Figure captions 

Table Baryon density (in fm~ 3 ), coupling constants, a and u> meson masses 
(in MeV) in our calculation. 

Fig.l Feynman diagrams for the calculation of scattering problem based on 
BNSD method. Double solid (dotted) lines represent exact nucleon(meson) 
propagators, and single solid(dotted) lines represent free ones. 

Fig.2 Feynman diagrams for the calculation of the self-energy. (a)Hartree 
diagram, (b)Fock diagram, (c) Higher order diagram more than fourth order. 
Doubly downward lines denote the hole states with momenta under the Fermi 
momentum kp and the doubly upward lines denote the intermediate states with 
momenta < k. 

Fig.3 The real part of nucleon self-energy extracted from Feynman propagator 
of nucleon on-shell, (a) The component of scalar part and (b) the component 
of vector part. 

Fig.4 The imaginary part of nucleon self-energy extracted from Feynman prop¬ 
agator of nucleon on-shell, (a) The component of vector part and (b) the com¬ 
ponent of scalar part. 

Fig.5 Optical potential in the Schrodinger equivalent form, (a) The imaginary 
part and (b) the real part. 

Fig.6 The nucleon mean free path in nuclear matter. The solid dots with error 
bars are experimental results from Ref. [1]. 

Table 


Pb 

m s 

m v 

9s 

9v 

0.170 

550 

783 

9.59 

11.67 

0.192 

550 

783 

9.12 

11.00 
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